Producing Bose condensates using optical lattices 
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We relate the entropies of ensembles of atoms in optical lattices to atoms in simple traps. We 
then determine which ensembles of lattice-bound atoms will adiabatically transform into a Bose 
condensate. This shows a feasible approach to Bose condensation without evaporative cooling. 



(N ■ 

o : 

o . 



(N 



> 

o 

O 



o 



There is considerable interest in combining weakly in- 
teracting Bose Einstein condensates (BECs) [Q with op- 
tical lattices |2j . The BEC-lattice combination is a dream 
model system: interactions among atoms in the BEC can 
be simply parameterized, and the optical lattice potential 
is exactly calculable. Work to date, both experimental 
and theoretical, has started with a BEC, and then put 
it in a lattice. For example, BECs in ID lattices have 
been used to make a pulsed atom laser Q and to demon- 
strate four wave mixing |^ and squeezing of atoms. 
The superfluid-Mott insulator transition has been pre- 
dicted § and observed § with BEC in a 3D lattice. In 
this paper, we theoretically address the converse prob- 
lem. Atoms that have never been Bose condensed start 
out deeply bound in a 3D optical lattice. The optical lat- 
tice is adiabatically removed, so that the atoms are left in 
either a flat-bottom or a harmonic trap. We determine, 
via entropy comparison, those initial conditions for which 
the final state is a BEC. Apart from the insight this pro- 
vides for the types of experiments performed to date, it 
also presents a strategy for achieving BEC without evap- 
orative cooling. This novel approach to BEC could be 
orders of magnitude faster than existing techniques. 

We consider only initial lattice distributions with ei- 
ther one or zero atoms per site. These are the most likely 
distributions after laser cooling in a lattice, because pho- 
toassociative collisions efficiently remove pairs of atoms 
at the same site We assume that atoms have been 
optically pumped and are thermally distributed among 
lattice site vibrational levels, the result of, for instance, 
3D Raman sideband cooling in an optical lattice |^ [0 . 

To realize these initial lattice distributions, we assume 
that atoms are initially trapped with no coupling between 
the sites and an infinitely strong in-site repulsion. Each 
site is represented by a three-dimensional harmonic os- 
cillator of frequency w. We also assume a site-dependent 
energy offset ej — ^hui of the bottom of each harmonic 
oscillator 



Hi = iei + hui ja\ fii, 

i {j} \ a=x,y,z / 



{,}+e(ni) , (1) 



where the operator = corresponds to 

the number of atoms in the i'th site in the "in-site" state 



given by the {j} = {jx, jy, jz} set of the "Cartesian" 
quantum numbers, the operator nj = gives 
the total number of atoms in the i'th site, the function 



e(n) 



0, for n = 0, 1 
+CXD for n > 2 



(2) 



provides the infinite in-site repulsion, and i^l {j}) 

is an operator creating (annihilating) a bosonic atom in 
the i'th site in the "in-site" state 

In what follows we will use a Grand-Canonical ex- 
pression for the entropy: S ~ — rr[ln(p)p] , where 
p = Z^^ exp(— /3(iJ — iiN)) is the Grand-Canonical 
A^-body density matrix of the system normalized to 
unity,Z ~ Tr[exp(— /3(7J — ^))] is the partition function, 
the chemical potential /j, is chosen to provide the actual 
total number of particles N: Tr[N p] = N , N being the 
total- number-of-particles operator, and /3 = 1/T is the 
inverse temperature corresponding to the total energy E 
of the system: Tr[H p] = E. 

Generic lattice: entropy. The total entropy of a lattice 
is the sum of the entropies of individual sites: 



■At 



(3) 



that according to expression for the entropy given above 
can be shown to be 



Si = In(zi) + l3{eh.o.{d) ~ pijn^ . 



(4) 



Here /ii = /i — Ci is the effective chemical potential, 
Zj = 1 + e^^^Zh^oXd) is the site's partition function, and 
Zh.o.id) = (1-e-^''")-'^ and eh.o.{d) ^ dhuj (e'^'"^-!)"! 
are respectively, the partition function and energy of a 
single d-dimensional harmonic oscillator. 

For a given temperature and chemical potential the 
total number of particles in the lattice is A^ = ni , 
where 



rii 



l^eP^'■' Zh.o.id) 
is the mean occupation of the site. 



(5) 
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Flat lattice: entropy. Here we will assume that all the 
lattice sites have the same energy offset and the number 
of available sites is iV, : 



ei = Yie[-m/2;+m/2f 



(6) 



The computation of the entropy is lengthy but straight- 
forward. It turns out that the entropy per particle 
si = Si/N (where here and below the index I stands 
for "lattice") is a universal function of the mean lattice 
site occupation, n = (the same for each site), and 
the probability, wo, that the ground vibrational state of 
a given site is occupied if the site itself is occupied (also 
the same for each site). The entropy-per-particle is 



si{wo,n) 



„T=0 



("-) + Sh.oXwo,3) 



(7) 



where 



sf="(n)=n-i [(l-n)ln(l/(l-n))+nln(l/n)] (8) 

is the entropy-per-particle in the same lattice but at zero 
temperature [ pT| , and 



i/d 



(9) 



wl^" Hl/wl^') + (1 - wl^') ln(l/(l - wl^")) 



is the entropy of a single one-particle d-dimensional har- 
monic oscillator whose ground state is occupied with a 
probability wq. Thermal equilibrium is assumed. 

Lattice with a harmonic envelope: entropy. Consider 
now a lattice whose site ground state energies harmoni- 
cally depend on the site position: 



(10) 



where I is the distance between the sites assumed to form 
a cubic lattice. The computation of the harmonic lattice 
entropy is pretty analogous to the one for the flat lattice. 
It leads to 



si{wo,n) ^ sf °{n) + Sh.o.{wo,^) , 



(11) 



where n — Tr[ni=o p] now denotes the site occupation in 
the center of the lattice, 



T=o^-^ 5Li5/2(-g(n)) 



2Li3/2(-9('^)) 



q{n) 



1 -n 



(12) 
(13) 



is the entropy-per-particle at zero temperature, and the 
d-dimensional harmonic oscillator entropy Sh.o.iwo, d) is 
given above. Here Li^(x) = ^"'/j'' is the polyloga- 

rithmic function. In the course of the above calculation 
we implicitly assumed that the site occupation was a slow 



function of the site index, and thus the summation over 
the sites was replaced by an integral. 

A general power-law trap: entropy at Tc. At the end 
of the thought transfer process atoms are confined in a 
smooth trap. In what follows, we will assume that many 
trap levels are populated, so that the semiclassical ap- 
proximation applies. Within the semiclassical approxi- 
mation, the number of (non-condensed) particles and the 
total entropy are given by 



St^ j d^rd3p(7(r,p) 
respectively, where the expressions 



(14) 
(15) 



w 



(r,p) = {2tiK) 3^np(n|r,p) 
1 



g/3(/it(r,p)-p) _ I 



(16) 



and 



CT(r, p) = {2'Kn) ln(p(n|r, p)) p(n|r, p) 

/3(/it (r,p) - /i) 



ri=0 



ln(z(r,p)) 



(17) 



g/3(/it(r,p)-;i) _ I 

are the particle- and entropy phase-space densities, and 
p(n|r, p) = z(r, p)-ie-'^('''('''P)-^)" (18) 

OO 

z(r,p) = ^e-'3(''*('-'P)-^)" (19) 

ri=0 

are the distribution of the number of particles in a volume 
(27r?i)'^ phase-space cube centered at a point (r, p), and 
the corresponding partition function. As has been shown 
in ||l^ , Bose-Einstein condensation occurs as soon as the 
total number of particles exceeds the upper bound for 
the number of the non-condensed particles N* = N' 

given by its value at zero chemical potential jT^ : 



N > N* 



BEC 



(20) 



Below we will reformulate this condition in terms of the 
entropy per particle St — St/N . 

Consider a general power-law trapping potential 



U{v)= «"(^")' 



(21) 



a—x.y.z 



After lengthy but straightforward calculations one can 
show that in this case the critical entropy per particle 
assumes a universal form 



N* 



(r; + 2)C(7? 2) 
C{V + 1) 



(22) 
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where 77 = 1/2 + J2a=x.y,zi''^a)~^ and ((x) is Rieman 
zeta-function. For the cases of a flat-bottom box and 
three-dimensional harmonic osciUator, rj is equal to 1/2 
and 2 respectively, leading to 



3D box: 



3D HO: 



5C(5/2) 
2C(3/2) 



1.284. 



= ^ = 3.602.. 



(23) 



Finally, assuming the trap has a positive fixed- heat 



capacity. Cm = -§f 



E 



N 



^ dT 



S > Q and the critical 



N 



number of particles monotonically increases with temper- 
N* > 0, one can show that the BEG condition 



ature 



d 

' dT 



( pO| ) is equivalent to the following condition on the en- 
tropy per particle: 



St < s* 



BEG 



(24) 



This form of the BEG criterion is particularly useful for 
analysis of adiabatic processes. 

Lattice-to-trap transfer process. Gonsider the following 
process. First, atoms are loaded into the optical lattice 
and thermalized. Next an external trap is built around 
the lattice. Finally, the original lattice is removed, and 
the state of the resulting gas is detected. The process is 
schematically depicted at Fig.|[ 



where s;(n, wq) is the entropy per particle in the original 
lattice, and is the critical entropy per particle for the 
given trap. 

The gas in the final trap will (not) be a BEG if the lat- 
tice entropy is below (above) the critical value s^. The 
lattice and critical trap entropies were explicitly com- 
puted above for two particular cases each. The four re- 
sulting phase diagrams are shown in Figs. and Bl 
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FIG. 1. An artist's view of the lattice-to-trap adiabatic 
transfer process. 

The transfer from lattice to trap is assumed to be ther- 
modynamically adiabatic, and therefore the entropy of 
the system is conserved throughout the process. In par- 
ticular, the entropies of the final trap and original lattice 
must be the same: St = Si. Since the number of particles 
N is also conserved, the above conservation law can be 
reformulated as 



St = si 



(25) 



where s = S/N is the entropy per particle. Gonsider 
a two-dimensional plane whose coordinates are given by 
the peak lattice site occupation n and ground state frac- 
tion in the occupied sites wq. The no-BEG vs. BEG 
transition line in this plane will be thus given by an im- 
plicit equation 
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FIG. 2. Phase diagram for flat-lattice-to-trap adiabatic 
transfer. The final traps are (a) flat-bottom box and (b) 3D 
spherically symmetric harmonic oscillator. The experimental 
points correspond to results in references [8] (circle) and [10] 
(diamonds). 

In a harmonic trap, the peak phase space density ex- 
ceeds the average phase space density. Because the peak 



3 



phase space density determines the BEC threshold p4[ , 
BEC is more readily obtained when lattice atoms are 
transferred to a harmonic trap than to a flat-bottom 
trap (compare Figures 2B and 3B to 2 A and 3A). When 
the initial distribution in the lattice is flat (Figures 2A 
and 2B), the initial bulk phase space density is constant 
across the sample. There is then more to be gained by 
changing the bulk shape of the trap than when the initial 
phase space density is already non-uniform (Figures 3A 
and 3B). The combination of these effects qualitatively 
explain why the flat lattice to harmonic trap scenario is 
the most promising approach to BEC (Figure 2B). 
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FIG. 3. Phase diagram for harmonic-lattice-to-trap adi- 
abatic transfer. The final traps are (a) fiat-bottom box and 
(b) 3D spherically symmetric harmonic oscillator. 



Fortunately, the photoassociative collisions that limit 
the average site occupation to one half also flatten the 
large scale distribution Q . Existing techniques can prob- 
ably produce transient site occupations of several, which 
will naturally lead to approximately flat top distributions 
after laser cooling. Some published experimental results 
for laser cooling in 3D optical lattices are marked on 
Figure 2B. Further increases in the ground state occu- 
pation are feasible especially by scattering light in 
the regime where the heating from rescattered photons is 
minimized [|l5|. As illustrated in Figure 2B, only a mod- 
est improvement is required to pass the BEC threshold. 

In summary, we have related the entropy of atoms dis- 
tributed in optical lattices to the entropy of atoms in 
harmonic and flat-bottom traps. By entropy comparison, 
we have calculated BEC phase diagrams for lattice-bound 
atoms. These diagrams create a framework for consider- 
ing experiments with BECs in optical lattices and suggest 
a way to obtain BEC without evaporative cooling. 
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